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Abstract 

We briefly report about a possible settlement of the still ongoing dispute concerning 
the existence of SUSY signals in 4jet events at LEP1. We base our arguments 
on a simple selection strategy exploiting secondary vertex tagging and kinematical 
constraints, which could allow one to access or exclude gluino events for a broad 
range of masses and lifetimes. 
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Introduction 



If the results of the LEP1 measurements were pieces of a jigsaw puzzle reproducing the 
edifice of the Standard Model, then one could well question that some of these are ap- 
parently not so perfectly shaped, i.e., they fit into their original location only with some 
effort. This is certainly the case for the determination of T c and r&, the partial widths of 
the Z into c and b quarks, for which claims of Supersymmetry (SUSY) hints have been 
made |]J. Indeed, there is another controversy still open along the same lines. Somewhat 
less glamorous but not for this less important is the possibility of gluino events being 
present in 4-jet decays [||. [|. 

The story goes as follows. The colour factors Ca, Cf and Tp of QCD || can be 
measured by fitting some angular distributions^ whose shape significantly depends on the 
partonic composition of the a 2 s decays of the Z. Then one obtains that, although the 
experimental measurements agree well with ordinary QCD, it is not possible to rule out 
its Supersymmetric version, which predicts that light gluinos g (the SUSY partners of 
the gluons) can be produced at LEP1 energies 0. In detail, gluinos with a mass rrig ~ 2 
GeV yield an expectation value for Tp/Cp that is within one standard deviation of the 
measured one ||. In fact, if the gluino mass is light enough ||, such particles should 
be produced in the process e + e~ — > QQgg via a g* —> gg splitting [0, ||. Since gluinos 
are coloured fermions, such events would enter into the sample with a behaviour similar 
to that of QQqq events. Naively, one could well say that the total number of flavours 
Np of the theory is apparently increased, such that, a SUSY signal reveals itself as an 
enhancement of Tr = NpTp. Such kind of new particles are at present still compatible 
with the experiments [^, [10[| . 

The reason why experiments have not given a conclusive answer so far is that both 
systematical (hadronisation, higher order perturbative corrections) and statistical (4-jet 
decays constitute only ~ 10% the hadronic sample) errors spoil considerably the precision 
of the measurements, thus preventing one from putting stringent bounds on Ca, Cf and 
Tr. However, the most serious and intrinsic limitation of the analyses performed up to 
now is that they made use of energy ordering to distinguish between quark and gluon jets 
and to assign the momenta to the final state partons^. 

A clear improvement to this approach is the one proposed in Ref. \T2\. There, it was 
shown the superiority of using 4-jet samples in which two jets are tagged as 6-jets. In 
this way, one gets a greater discrimination power between QQgg and QQqq events. First, 
because this way one is able to distinguish between quarks and gluons, thus assigning the 
momenta correctly. Second, because gluon rates are reduced by almost a factor of 2 with 
respect to the quark ones, such that differences between the two partonic components can 



1 That is, the angles of Bengtsson-Zerwas, of Korner-Schierholz-Willrodt, of Nachtmann-Reiter and 
that between the two least energetic jets 

2 The two most energetic jets are identified as primary quarks. Unfortunately, for QQgg events, in 
only half of the cases the two lowest energy partons are both gluons |0| ! 



be more easily investigated. Following Ref. |]T2|, the LEP Collaborations have recently 
performed new studies Hl3| , [14}] , whose preliminary results show indeed that SUSY predic- 
tions can be more efficiently constrained. Furthermore, they have proved that adopting 
a double vertex tagQ does not ruin the advantages gained with particle identification. 

Besides the final results of these new, improved analyses, we want to stress that there 
are other possibilities offered by the //-vertex devices, that can be exploited in order to 
either confirm or disprove the presence of SUSY signals in 4-jet events. This is apparent 
if one notices that light gluinos can also be relatively long-lived, such that they might 
produce detectable secondary vertices j0|. It is the purpose of this letter to study to 
which extent such experimental techniques can be used for detecting or ruling out SUSY 
signals at LEP1, even when no special effort is made to distinguish between displaced 
vertices due to 6-quarks and to gluinos. 

The plan of the paper is as follows. In the next Section we describe our calculations, 
in Section 3 we discuss the results, and in Section 4 we summarise and conclude. 

2. Calculation 

In carrying out the study described here we made use of the FORTRAN matrix elements 
already discussed in Ref. [0 and presently used for experimental simulations 0, upgraded 



with the inclusion of the gluino production and decay mechanisms (see also Ref. [17|)n- 
The programs do not contain any approximations, the intermediate states 7* and Z being 
both inserted, and the masses and polarisations of all particles in the final states (of the 
two-to-four body processes) retained. The availability of the last two options is especially 
important if one considers, on the one hand, that in 6-tagged samples all final states are 
massive, and, on the other hand, that in proceeding to experimental fits one could well 
select restricted regions of the differential spectra of the angular variables, where the rates 
are likely to strongly depend on the spin state of the partonsf]. 



As jet finding algorithm we have adopted here the Durham (D) scheme []18| |. However, 
none of the results drastically depends on the choice of the jet recombination procedure 
and/or the value of the jet resolution parameter, y cut . Finally, to make clear the rest of 
the paper, we use the following notation: when heavy flavour identification is implied, 
labels 1 & 2 refer to the two tagged jets and 3 & 4 to the two remaining ones. If no vertex 
tagging is assumed, jets are labelled according to their energy, E\ > E 2 > E 3 > E4. 



3 Which reduces considerably the statistical sample, as the current efficiency at LEP1 in tagging a 
displaced vertex is e ~ 30% per jet. 

4 The numerical values adopted for quark masses and SM. parameters can be found in Ref. JH| . 

5 For this reason we have not used the results published in literature for the gluino decay rates, as 
these are averaged over the helicities of the unstable particle. Instead we have recomputed the relevant 
Feynman decay amplitudes by preserving the gluino polarisation and by matching the latter with the 
corresponding one in the production process. 



3. Results 



3.1 Gluino tagging 

When dealing with tagging a secondary vertex possibly due to gluino decays, several 
points must be addressed. First of all, one has to confine oneself to secondary vertex 
analyses onlyQ, however this technique has a larger efficiency than any other method |B| . 
Second, the vertex has to be inside the detectors, so that only gluinos with r g < 10~ 9 s 
can be searched for ||T5| . Nonetheless, this represents an appealing opportunity, as a 
substantial part of the (m g , Tg) window [0] not yet excluded by the experimental data 
could be covered. The latest constraints still allow for the existence of relatively long- 
lived and light gluinos, in the parameter regions: (i) m g ~ 1.5 GeV and r g ~ 10~ 8 s; (ii) 
m g > 4 GeV and r g > 1(T 10 s |, [H|. 

In this respect, we exploit a sort of 'degeneracy' in lifetime between 6-quarks and 
gluinos, assuming that when making secondary vertex tagging one naturally includes in 
the 262jet sample also SUSY events, in which a g behaves as a b. We call such approach 
'minimal trigger' procedure, as we propose a tagging strategy that does not take into 
account any of the possible differences between gluinos and 6-quarks in 4-jet events with 
two secondary vertices (thus, in the following, we will generally speak of 'vertex tagging'). 
There are in fact at least three obvious dissimilarities. 

(i) Their charge is different, such that one could ask that the jet showing a displaced 
vertex has a null charge. This would allow one to isolate a sample of pure SUSY events. 
However, we remind the reader that measuring the charge of a low energy jet in 4-jet 
events would have very low efficiency (in isolating a very broad hadronic system in an 
environment with high hadronic multiplicity) and has not has not been attempted before. 

(ii) Gluino lifetimes much longer than ^-lifetimes are still consistent with experiment (note 
that Tf, R3 1CU 12 s), such that recognising a displaced vertex with decay length d 3> 3 mm 
(that of the b) would allow one to immediately identify gluinos. Unfortunately, most of 
the 6-tagged hadronic sample at LEP1 has been collected via a bi-dimensional tagging 
(see, e.g., Ref. PD[ ). Thus, different d's could well appear the same on the event plane. 
Furthermore, tagging a d > 3 mm vertex would allow one to separate gluinos with r g > 75, 
but this would not be helpful if r g < r b . 

(iii) Other than in lifetime, 6-quarks and gluinos can differ in mass as well, such that 
one might attempt to exploit mass constraints to separate SUSY and pure QCD events. 
However, on the one hand, one could face a region of mi,-m g degeneracy and, on the other 
hand, one should cope with the ambiguities related to the concept of mass as defined at 
partonic level and as measured at hadron level. 

We emphasise that measuring the charge of the vertex tagged jet, attempting to disentan- 

6 Thus neglecting other forms of heavy flavour tagging: such as the high px lepton method, as gluinos 
do not decay semileptonically. 



gle different decay lengths or measuring partonic masses could well be further refinements 
of the procedure we are proposing. These could be implemented at a later stage without 
spoiling the validity of our approach. In addition, all these aspects would necessaril need 
a proper experimental analysis, which is beyond the scope of a theoretical study. 

The steps of our procedure are very simple. Under the assumption that fe's and g's 
are not distinguishable by vertex tagging, one naturally retains in the sample all SUSY 
events, whereas the ordinary QCD components are reduced by a factor of 5 (QQgg) and 
3 (QQqq). Then, it is easy to notice that there exist clear kinematic differences between 
the QQgg, QQqq and QQgg components. This is shown in Fig. 1, where we plot the 
quantities = Mfj/s = (pi +pj) 2 /s, where ij = 12 or 34 and s = Mf. The behaviour 
of the curves is dictated by the fact that gluinos are always secondary products in 4- 
jet events, whereas quarks and gluons are not (lower plots). When no vertex tagging 
is exploited and the common energy ordering is performed, such differences are washed 
away (upper plots). The value m y = 5 GeV is assumed for reference, the shape of the 
distributions being qualitatively the same regardless of it. Therefore, if one simply asks to 
reject events for which, e.g., Y 12 > 0.2 and/or Y34 < 0.1, one gets for the total rates of the 
three components the pattern depicted in Fig. 2. Notice that the drastic predominance 
of 2Q2g events in the complete 'unflavoured' sample has disappeared. Further, the total 
rates of ordinary QCD events are significantly reduced compared to those of SUSY events 
(after vertex tagging, top right), and eventually the QQgg fraction is always comparable 
to that of QQgg + QQqq events (when also the kinematics cut are implemented, bottom 
left). Most important, this is true independently of the gluino mass, of the jet algorithm 
and of the y cut value used in the analysis. 

Therefore, after our event selection, SUSY signals would certainly be identifiable, as 
a clear excess in the total number of 4-jet events with two displaced vertices. Thus, the 
presence of gluinos at LEP1 could be revealed or excluded at least over appropriate regions 
of masses and lifetimes]]. We finally stress that, as we are concerning here with total rates 
and not with differential distributions, the event number should be sufficient to render 
the analysis statistically significant!. Furthermore, the ambiguities related to the fact 
that gluino effects on the total number of 4-jet events are comparable in percentage to 
the systematic uncertainties due to the jet hadronisation process and/or the y cut selection 
procedure are much less severe in ours than in the usual approach 0. However, since the 
key point of the present study is to exploit the b/g vertex degeneracy, a highly enriched 

7 And this should certainly be done after the appropriate MC simulations, including the details of the 
detectors and of the tagging procedure as well as a generator where trig and r g enter as free parameters 
to be determined by a fit. 

8 In this respect we acknowledge that many of the aspects of our approach were already employed in 
Rcf . jl5) , however the tagging procedure sketched there is well beyond the statistical possibilities of the 
LEP1 experiments. 

9 These are in fact the underlying difficulties of any analysis based on the 'unflavoured' hadronic sample 
and/or the jet energy ordering, which have not been overcome even in recent improved analyses [ pi) . 



heavy flavour sample should be selected in this case. 
3.2 Gluino decays 

Before closing, we should mention that a further aspect must be kept into account when 
attempting our analysis. It concerns the kinematics of the gluino decays. In the most 
widely supported SUSY framework [P2j] , the dominant gluino decay modes are g — > qq'y 
and g —>■ g^y, where 7 represents a 'photino' (better, the Lightest Supersymmetric Particle, 
which is a superposition of the SUSY partners of the neutral gauge bosons of the theory). 
Furthermore, the qq'y channel is, in general, largely dominant over the #7 mode ||23|| . 

The crucial point is that in both cases the gluino decays into a jet with missing energy. 
It is not our intention to discuss the possibility of selecting such a signature, as we are 
mainly concerned here with the fact that the energy left to the hadronic system Eh is 
above the experimental cuts in minimal hadronic energy, which are used to reduce the 
backgrounds (e.g., in Ref. || the threshold was set equal to 3 GeV). In Fig. 3 (first 
three plots) we show the Eh spectra after the gluino decay, in both the channels. Two 
kinematic decay configurations are considered: a massless photino, and a massive one 
(i.e., = l/2mg). The message is that in the most likely SUSY scenario (i.e., three- 
body decay dominant and massless photino) all gluino events should be retained in the 
event selection. Conversely, the figure illustrates the percentage of these which will pass 
the adopted trigger requirements. Finally, in the bottom right plot of Fig. 3 we show 
the dependence of the SUSY rates on the value of m- g . Below rrig m 5 GeV, the mass 
suppression is always less than a factor of 2. 

4. Summary and conclusions 

In this paper we have stressed the importance of using at LEP1 samples of 4-jet events 
in which two of the jets show a displaced vertex, and of adopting simple invariant mass 
cuts based on the different kinematics of partons in the final state. This can help to settle 
down the ongoing dispute about the existence of SUSY events in the data, at least for a 
wide range of gluino masses and lifetimes. Those presented are theoretical results, which 
should be in the end verified by detailed MC simulations that could even improve our 
event selection strategy (by exploiting differences between 6's and <fs, in charge, mass 
and lifetime). Therefore, it is our opinion that the matter raised and procedures similar 
to the ones outlined here would deserve experimental attention. An enlarged and more 
detailed version of the present paper, which contains a generalisation of our results to 
other three jet schemes together with a discussion of the angular variable dependence of 
the ordinary QCD and SUSY rates, will be given elsewhere ||24|| . 
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Distributions in the rescaled invariant masses = M^/s, where ij = 12,34, for ordinary QCD 
and for SUSY 4-jet events, in the D scheme with y® ut = 0.002, without (Q = q) and with (Q = b) 
vertex tagging. Here, nig = 5 GeV. 
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Fig. 2 



Total cross sections of ordinary QCD and of SUSY 4-jet events, in the D scheme, without (Q = q) 
and with (Q = b) vertex tagging, and after the kinematic cuts, for three different values of rrig. 
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Differential distributions in the hadronic energy of the 'gluino jet' after the two possible SUSY 
decays, in the D scheme, for various combinations of m g and m^; and total cross section of 
gluino events in 4-jets, as a function of m g and for three different values of y^ut- 



